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Shock  tunnel  operation  and  correlation  of 
boundary  layer  transition  on  a  cone  in 
hypervelocity  flow 


J.S.  Jewell1,  J.E.  Shepherd1,  and  I. A.  Leyva2 


1  Introduction 

The  Caltech  T5  reflected  shock  tunnel  is  used  to  produce  hypervelocity  flow  over  a 
range  of  velocities  and  pressures  by  varying  the  test  gas  and  operating  parameters 
of  reservoir  enthalpy  (hres)  and  reservoir  pressure  ( Pres )•  One  area  of  research  in  T5 
is  the  measurement  of  boundary  layer  behavior  and  transition  from  laminar  to  tur¬ 
bulent  flow  on  a  smooth  5 -degree  half-angle  cone  [3,  1,  11].  To  design  experiments 
that  involve  the  measurement  or  manipulation  of  instability  and  transition  processes 
(for  example,  Jewell  et  al.  [7]),  it  is  important  to  choose  tunnel  conditions  for  which 
the  expected  transition  location  is  at  least  approximately  known.  In  the  present  pa¬ 
per,  we  discuss  the  selection  of  tunnel  operating  parameters,  the  correlation  of  those 
parameters  with  measurements  of  boundary  layer  transition,  and  some  observations 
on  the  analysis  of  transition  location  in  terms  of  local  boundary  layer  properties. 


2  Tunnel  Operation 


Flow  conditions  in  T5  are  calculated  from  three  tunnel  measurements:  the  shock 
speed,  initial  shock  tube  fill  pressure  and  composition,  and  reservoir  pressure  at  the 
end  of  the  shock  tube  during  the  run  time  [4].  Only  experiments  with  measured 
shock  speeds  that  fall  within  the  uncertainty  for  the  adjusted  shock  speed  curve 
predicted  by  the  shock  jump  conditions  from  the  primary  diaphragm  burst  pressure, 
driver  gas  composition,  and  initial  shock  tube  conditions  are  included  in  the  present 
data  set. 

There  are  a  number  of  other  potential  sources  of  measurement  error,  bias,  or  un¬ 
certainties.  These  include:  nonideal  gas  behavior  in  the  reservoir  due  to  the  high 
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pressure;  the  extrapolation  of  the  shock  speed  (which  decays  as  it  propagates  down 
the  shock  tube)  to  the  end  wall;  nonuniformity  of  reservoir  conditions  due  to  non¬ 
ideal  shock  reflection;  and  the  method  of  correcting  flow  conditions  from  the  ideal 
reflected- shock  pressure  to  measured  reservoir  pressure  using  an  isentropic  expan¬ 
sion.  Furthermore,  the  1-D  nozzle  computation  does  not  account  for  boundary  layer 
growth  within  the  nozzle,  off-design  operation  conditions  that  lead  to  flow  nonuni¬ 
formity,  or  vibration- translation  nonequilibrium  and  freezing  within  the  nozzle, 
which  is  significant  for  the  N2  cases.  For  the  uncertainties  that  can  be  quantified, 
we  have  combined  these  to  obtain  the  error  bounds  on  measured  properties  that  are 
shown  in  this  paper.  The  tunnel  parameters  for  the  present  studies  in  air  and  nitrogen 


Fig.  1  Tunnel  operating  parameters  hres  and  Pres  for  the  present  studies  in  air  and  N2,  compared 
with  past  results  from  Adam  [1]  and  Jewell  et  al.  [5]. 


are  compared  with  those  of  two  past  data  sets  in  air  in  Figure  1 .  The  present  work 
both  overlaps  and  extends  the  parameters  of  the  past  studies,  especially  for  low  pres¬ 
sure  and  enthalpy.  The  R2  values  for  the  correlation  between  the  two  parameters  are 
respectively  0.52  and  0.65  for  the  present  N  2  and  Air  data  sets,  and  respectively  0. 10 
and  0.59  for  the  Adam  [1]  and  Jewell  et  al.  [5]  data  sets. 

Freestream  conditions  are  taken  as  the  conditions  at  the  nozzle  exit.  The  100:1 
area  ratio  contoured  nozzle  is  designed  to  operate  at  Mach  6.  Because  the  shape 
is  optimized  for  a  single  condition,  there  is  significant  variation  of  the  exit  Mach 
number  over  the  range  of  tunnel  operating  parameters,  presented  in  Figure  2  for  air 
over  the  conditions  of  the  present  study. 


3  Experiments 


The  experimental  model  is  a  1-m  long,  smooth,  5 -degree  half-angle  cone  with  a 
nominally  sharp  tip  of  radius  0.18  mm  and  oriented  at  a  zero  angle  of  attack.  The 
range  of  Reynolds  numbers  evaluated  at  the  boundary  layer  edge  and  Dorrance  [2] 
reference  temperature,  which  is  used  as  representative  of  conditions  within  the 
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Fig.  2  Calculated  nozzle  exit  velocity,  sound  speed,  and  Mach  number  in  air  over  a  range  of  tunnel 
operating  parameters  hres  and  Pres . 


Fig.  3  Reynolds  number  evaluated  at  the  boundary  layer  edge  in  air  (left)  and  at  Dorrance  reference 
conditions  in  air  (right)  over  a  range  of  tunnel  operating  parameters  hres  and  Pres . 


boundary  layer,  is  presented  for  air  in  Figure  3.  The  cone  is  instrumented  with  80 
thermocouples,  providing  heat  transfer  measurements  from  which  transition  loca¬ 
tion  may  be  determined.  Heat  transfer  results  are  normalized  by  Stanton  number 
and  Reynolds  number,  and  the  location  of  transition  onset  determined  as  described 
in  Jewell  et  al.  [7]. 

Parameters  for  experiments  that  bracket  the  range  of  conditions  studied  in  both 
air  and  N2  are  presented  in  Table  1.  A  supplemental  report  [6]  is  available  online, 
which  includes  the  full  tables  and  plots  for  all  34  experiments  included  in  the  present 
study. 


4  Results  and  Analysis 

Results  for  the  location  Xtr  of  transition  onset  on  the  cone  are  presented  in  terms  of 
hres  and  Pres  in  Figure  4  for  both  N2  and  air.  Multivariable  linear  regression  analysis 
is  performed  on  these  data  sets  (for  details,  see  [6]).  Both  the  present  N2  and  air 
results  have  a  positive  dependence  on  hres  (linear  model  coefficient  of  0.56  for  N2, 
0.55  for  air)  and  a  negative  dependence  on  Pres  (linear  model  coefficient  of  —0.45 
for  N2,  —0. 15  for  air).  The  historical  air  data  of  Adam  and  Hornung[l],  are  analyzed 
in  the  same  way,  and  likewise  show  a  significant  positive  dependence  of  Xtr  on  hres 
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Table  1  Parameters  for  selected  experiments  spanning  the  present  range  of  conditions  for  air  and 
N2. 


Gas 

hres 

[MJ/kg] 

p 

1  res 

[MPa] 

Re/m 

[1/m] 

Re*  /m 

[1/m] 

Xtr 

[m] 

$99tr 

[mm] 

f,r 

[kHz] 

2744 

Air 

7.68 

60.7 

8.19  x  106 

5.39  x  106 

0.51 

0.98 

1094 

2758 

Air 

11.07 

72.0 

6.11  x  106 

4.79  x  106 

0.74 

1.26 

1002 

2764 

Air 

5.27 

16.5 

3.62  x  106 

1.83  x  106 

0.52 

1.65 

551 

2773 

n2 

8.99 

16.7 

2.02  x  106 

1.25  x  106 

0.67 

2.26 

522 

2776 

n2 

7.17 

45.9 

7.09  x  106 

3.94  x  106 

0.39 

0.97 

1102 

2783 

n2 

15.88 

53.3 

3.20  x  106 

2.62  x  106 

0.63 

1.56 

966 

(linear  model  coefficient  of  0.72)  and  negative  dependence  on  Pres  (linear  model 
coefficient  of  —0.28). 

Both  the  present  N2  and  air  results  have  a  positive  dependence  on  Pres  (linear 
model  coefficients  of  0.31  for  N2,  0.59  for  air)  for  the  transition  Reynolds  number 
evaluated  at  Dorrance  reference  conditions,  Re  *r,  but  neither  have  a  dependence  on 
hres  that  is  statistically  significant.  The  historical  air  data  of  Adam  and  Homung[l] 
likewise  show  a  significant  positive  dependence  of  Re  *r  on  Pres  (linear  model  coeffi¬ 
cient  of  0.34),  but  no  statistically  significant  dependence  on  hres. 

The  usual  approach  [12]  for  representing  boundary  layer  transition  onset  is  in 
terms  of  a  transition  location  Reynolds  number  Retr  or  Re*r  and  previous  analy¬ 
sis  [1]  of  T5  data  has  utilized  this  approach.  However,  in  hypervelocity  flow  with  a 
cold  wall,  the  principal  boundary  layer  instability  mechanism  is  the  predominantly 
inviscid  acoustic  or  Mack  mode  [9].  Unlike  the  viscous  instability  of  low-speed 
boundary  layers,  at  high  speeds  the  role  of  viscosity  is  primarily  in  determining  the 
mean  flow.  The  properties  of  the  acoustic  instability  are  determined  by  the  local 
boundary  layer  thickness  and  profiles  of  velocity  and  thermodynamic  properties. 
This  suggests  the  approach  of  correlating  transition  distance  with  Xtr/  5  99  and  the 
acoustic  properties  of  the  boundary  layer  rather  than  a  Reynolds  number. 

As  shown  in  Figure  5(1),  Xtr/  699  is  relatively  independent  of  edge  Mach  number 
but  shows  a  systematic  dependence  on  the  gas  type  and  pressure.  The  scaled  dis¬ 
tance  ranges  from  about  300  for  low  pressure  tests  up  to  600  for  high  pressure.  The 
radiation  of  acoustic  disturbances  from  the  turbulent  boundary  layer  on  the  nozzle 
wall  and  the  jet  shear  layer  is  an  important  source  of  disturbances  in  ground- testing 
facilities  [12],  and  varying  receptivity  of  the  boundary  layer  to  these  radiated  distur¬ 
bances  is  a  likely  explanation  for  the  trends  observed  in  Figure  5. 

Lowering  the  pressure  creates  a  thicker  boundary  layer,  and  therefore  lowers 
the  most  amplified  second  mode  frequency  /  «  0.6Ue/28,gg.  This  may  account 
for  the  striking  correlation  of  Xtr/  5  99  with  the  frequency  ftr  that  is  shown  in  Fig¬ 
ure  5(r).  The  influence  on  transition  location  can  be  explained  by  the  measurements 
of  Parziale  et  al.  [10],  who  showed  that  in  T5,  most  of  the  noise  in  the  free  stream  is 
at  relatively  low  frequencies  (<  500  kHz),  and  observed  a  decrease  in  rms  density 
fluctuations  with  increasing  frequency.  This  is  consistent  with  the  present  observa- 
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tions  of  earlier  transition,  at  lower  most-amplified  frequencies,  for  lower  values  of 
Pres  in  both  air  and  N2. 


Fig.  4  (1)  Transition  onset  location  Xtr  for  N2,  and  (r)  air  in  terms  of  reservoir  enthalpy  hres  and 
pressure  Pres. 


Fig.  5  (1)  Scaled  transition  distance  Xtr/8,ggtr  vs.  Medge,  the  boundary  layer  edge  Mach  number, 
(r)  Scaled  transition  distance  Xtr/S^tr  vs.  ftr,  the  approximate  most-amplified  second  mode  fre¬ 
quency  at  transition. 


5  Conclusion 

We  have  re-examined  the  correlation  of  transition  onset  data  with  tunnel  operat¬ 
ing  conditions  by  modeling  transition  onset  location  with  two-dimensional  response 
surfaces  in  terms  of  hres  and  Pres  variables.  We  observe  a  positive  correlation  of  Xtr 
with  hres  in  both  air  and  N2  and  a  negative  correlation  of  Xtr  with  Pres.  The  pa¬ 
rameter  Re*tr  exhibits  a  positive  correlation  with  reservoir  pressure  in  both  gases. 
Controlling  for  variations  in  Pres ,  no  statistically  significant  dependence  of  Re*r  on 
hres  was  found  for  either  the  present  air  or  N  2  data  or  the  historical  air  data  in  Adam 
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and  Homung[l].  We  explore  an  alternative  normalization  of  the  transition  onset  lo¬ 
cation  by  the  local  laminar  boundary  layer  thickness  Xtr/8m 99  and  find  that  this  is 
essentially  independent  of  the  edge  Mach  number  for  a  given  gas  type  and  pres¬ 
sure  range  (Figure  5).  We  examine  the  correlation  of  transition  onset  with  the  most 
amplified  acoustic  frequency  in  the  boundary  layer  and  find  that  lower  frequencies 
correlate  with  smaller  values  of  normalized  transition  distance  Xtr/ 699,  suggesting 
that  the  frequency-dependent  amplification  of  the  tunnel  noise  may  be  responsible 
for  the  observed  systematic  variations  in  transition  onset  distance. 
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